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ABSTRACT 

We report on the RR Lyrae variable star, MACHO 176.18833.411, located toward the Galactic bulge 
and observed within the data from the ongoing Bulge RR Lyrae Radial Velocity Assay (BRAVA-RR), 
which has the unusual radial velocity of —372 ± 8 km s“^ and true space velocity of —482 ± 22 km s“^ 
relative to the Galactic rest frame. Located less than 1 kpc from the Galactic center and toward a 
field at (1,6)=(3,—2.5), this pulsating star has properties suggesting it belongs to the bulge RR Lyrae 
star population yet a velocity indicating it is abnormal, at least with respect to bulge giants and red 
clump stars. We show that this star is most likely a halo interloper and therefore suggest that halo 
contamination is not insignificant when studying metal-poor stars found within the bulge area, even 
for stars within 1 kpc of the Galactic center. We discuss the possibility that MACHO 176.18833.411 
is on the extreme edge of the bulge RR Lyrae radial velocity distribution, and also consider a more 
exotic scenario in which it is a runaway star moving through the Galaxy. 

Subject headings: Galaxy: bulge; Galaxy: kinematics and dynamics; Galaxy: structure; surveys 


1. INTRODUCTION 

It is well known that within the Galactic bulge at lon¬ 
gitudes |Z|<10°, there is a bar-like structure with a bar 
angle in the range 20 - 30°, traced by numerous old 
stellar population probes (e.g., red clump giants, Stanek 
1994, Miras, Groenewegen & Blommaert 2005). This 
bar/bulge is a rotating Box/Peanut (B/P) structure with 
a X-shape protrusion, made up largely of old and metal- 
rich stars (~10 Gyr, [Fe/H] falling between —0.5 and 
-1-0.5 dex, e.g., Johnson et al. 2013; Gonzalez et al. 2012; 
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Wegg & Gerhard 2013). 

There is speculation that the Milky Way also has an 
older, more spheroidal bulge population, and perhaps 
the greatest possibility of uncovering such a component 
would be within the most metal-poor bulge stars. One 
example of a metal-poor bulge population is the RR 
Lyrae stars (RRLs), and over the past few years, ^38 000 
RRLs toward the bu lge have been identifie d from photo¬ 
metric surveys (e.g., Soszyhski et al.|20l4 ). These RRLs 
are thought to exhibit a small metal ficity spread and are 


centered around [Fe/H] = —1 de x (Walker &: Terndrup 


1991 Kunder & Ghaboyer 2008 Soszyhski et al. 2014). 


I'he small [I'e/HJ spread rnay suggest that RRLs trace 
a more ancient stellar population than the majority of 
bulge RGB and red clump stars, which are more metal- 
rich, on average, and likely do not evolve t o become RRLs 
(e.g. Lee||1992 Walker fc Terndru^|1991 ). 

Recent photometric studies have reached different con¬ 
clusions regarding the relationsh i p of t he bulge RRLs 
and the bar/bulge. Dekany et al. (2013) combined opti¬ 
cal and infrared photometry of ^8000 ss OGLE-Hl dis¬ 
covered RRLs to find that unlike the red clump giants 
(RCGs), RRLs do not trace a strong bar. Instead, they 
have a more spheroidal, centrally concentrated distribu¬ 
tion, indicating that the RRLs belong t o a classical bulge 


that has co-evolved with the bar (e.g., Saha & Gerhard 


2013). 


Tn contrast, using 28 000 RRLs fro m the more spatially 
extend ed OGLE-IV bulge sample, [Pietrukowicz et aL] 
(2014) assert that the RRLs trace closely the barred 
structure formed of RCGs, and hence that the bulge 
RRLs are in the same gravitational potential together 
with the more massive Galactic bar. 

Missing still are the kinematics of the bulge RRLs, 
which can resolve this discrepancy and provide an un¬ 
derstanding of the origin of the old, metal-poor bulge 
component. Th e last pub l ished paper on bulge RRLs 
radial velocities, Gratton (1987), used a sample of 17 
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RRLs to conclude that the kinematic properties of RRLs 
in Baade’s Window are similar to that of the Miras, M- 
giants, K-giants, OH/IR sources and planetary nebulae 
of the Galactic bulge. 

In this paper, we report on a high-velocity RRL found 
serendipitously in the Bulge RR Lyrae Radial Velocity 
Assay (BRAVA-RR) to have 1 ;,.=—372 km s“^, which 
is well above the typical speed of the stars one might 
expect to find in the bulge. High-velocity stars are in¬ 
triguing in part because they can provide insight to the 
mechanisms that produce their velocities. The origin of 
high-velocity stars can also provide useful information 
about the environments from which they are produced. 
Here we investigate the cause for the high-velocity of MA¬ 
CHO 176.18833.411 to discern whether it is consistent 
with stars in the Galactic bulge and what this suggests 
about the formation of the Galaxy. 


2. OBSERVATIONS AND RADIAL VELOCITY 

MACHO 176.18833.411 was originally catalogu ed by 
the MACHO survey as a fundamental mode RRL (Kun¬ 


der et al. 2008[ | 




The star was one of the ~100 
RRLs surveyed spectroscopically as part of BRAVA-RR 
(NOAO PropID: 2014A-0143; PI: A. Kunder) in a field 
at (/, 6 ) = (3,—2.5) using the AAOmega multifiber spec¬ 
trograph on the Anglo-Australian Telescope (AAT). We 
observed this star twice on June 21, 2014, separated in 
time by 7 hours, and the observations were taken in dual 
beam mode centered on 86 OOA, with the 580V and 1700D 
gratings to probe the Calcium Triplet. This covers the 
optical window from about 8300A to 88 OOA at a resolu¬ 
tion of R~10,000. 

The data were reduced using the automated pipeline 
supplied by AAOmega, 2DFDR, and the spectra were 
cross-correlated using the IRAF cross-correlation rou¬ 
tine, xcsoa. Four Bulge RAdi al Velocity Assay 
(BRAVA) stars (Kunder et al. 2012) observed with the 
same setup were selected as radial velocity standard 
stars. Due to less than optimal weather conditions, the 
signal-to-noise ratio is ~i0 (see spectrum in Figure!^, 
and the consistency of our velocity result is 8 km s^, 
in agreement with the errors reported by xcsao for each 
individual measurement. 

Figure^ (middle) shows the pulsation curve using the 
radial velocity te mpla t e and scaled by its V-amplitude 
as outlined from Liu (1991). The two radial velocity 
measurements were folded by the known period to find 
the radial velocity as a function of phase. Both mea¬ 
surements fit the radial velocity template well for a RRL 
with a line-of-sight radial velocity, Vios of —372 km s“^. 
It is noteworthy that we did not adjust the template in 
phase, which indicates that the OGLE time of maximum 
brightness is reliable for this star. 


3. THE PROPERTIES OF MACHO 176.18833.411 

The OGLE-IV catalogue of RRLs provides V- and I- 
band light curves of their stars, with the photometric 
observations spanning from March 2010 to October 2013 
(Soszyhski et al.||2014|). The optical light curve for MA- 
CHO 176.18833.411 from the OGLE-IV observations is 
shown in Figure(left), and its most important proper¬ 
ties are summarised in Table [H 


This RRL is also designated OGLE-BLG-RRLYR-10353. 


To find the distance to the RRL, we first use the mean- 
flux magnitude as listed by OGLE-IV. Next, the EfV-IJ 
color excess along the stars line of sight is calculated from 
the observed V-I color at minimum light, (V-I) 


ohs 


(see e.g., Guldenschuh et al. 2005), by carrying out a 
Fourier fit to the V and l-band light-curves. 

The extinction, Aj can then be derived using: 


Ai = 0.7465E(V - I) + 1.37E{J - K), 


( 1 ) 


as introduced in Nataf et al.| (2013). Here , E(J-K) 


was taken fr om the bulge reddening maps of Gonzalez 
et al. (2012). Finally, it has been shown that RRLs 


F-bahd light curves and the phase difference of their 
Fourier decomposition can be used to estim ate photo¬ 


metric metall icities to an accuracy of 0.2 dex( Jurcsik & 
Kovacs 1996), although with high-qualit y light curves 


the fitting accuracy can be ^0.12 dex (Kovacs 2005). 
This relation is well-calibrated over the metallicity range 
from 2.0 dex to ^0 dex. Applying this method to 
the OGLE V-band light curve, the metallicity ([Ee/H]) 
of MACHO 176.18833.411 is estimated and placed on the 
Carretta et al. (2009) metallicity scale. Using the recali- 
bratiq n of the RRL luminosity scale by|Catelan fc Cortes] 


(2008), MACHO 176.18833.411 has an absolute magni¬ 
tude of My=0.58±0.13 mag. Similarly, using a quadratic 
relation between RRLs absolute magnitude and metallic¬ 
ity from Bono, Caputo, & di Criscienzo (2007), we find 
My=0.61±0.08 mag, where 0.08 is a reasonable error in 
the RRLs absolute magnitude-[Ee/ H] zero-point calibra¬ 
tion. The Benedict et al. (2011) My-[Fe/H] relation, 
however, results in iWv=U.43±U.07, indicating that the 
level of agreement between independent My measure¬ 
ments is as large as ~0.2 mag. A systematic uncertainty 
of 0.2 mag My is factored into the uncertainty in the es¬ 
timated distance, and does not change the results signifi¬ 
cantly. The apparent magnitude, reddening and absolute 
magnitude of MACHO 176.18833.411 lead to a distance 
of (m — M)o=14.29±0.13 mag. 


In the OGLE-HI catalog, Soszynski et al. (2011) 
marked ^400 stars as possessing proper motions relative 
to mean motion of the bulge that are large enough to be 
easily detected. MACHO 176.18833.411 is included in 
this list. To determine its proper motion, the centroids 


of stars were measured on ea ch OGLE-IV (Udalski, Szy- 
manski & Szymanski ||2015 ) image separately. All mea¬ 
sured centroids were transiormed to the common grid 
using positions of bright red giants i.e., stars belonging 
to the bulge, not the disk population. The transformed 
positions were fitted with a model that takes into account 
proper motion as w ell as differential refraction effects (see 
Poleski et al.||2013 for details). 

With our proper motion, the space velocity and po¬ 
sition vector can be resolved, and the 6 D position and 
velocity information for MACHO 176.18833.411 is given 
in Table 1. Eor this convention, the Sun’s orbital ve¬ 
locity vector Vq= [[/ 0 ,V 0 , 1 F©] = [14.O,12.24,7.25] km s“^, 
1 A.sr= 220 km s“^ and position = [8.28,0,0] kpc. The de¬ 
rived true space velocity is —482 ± 22 km s“^ relative to 
the Galactic rest frame. 


4. DISGUSSION 

4.1. Possible explanations for the high-velocity star 














































BRAVA-RR 


3 



Phase 




Wavelength (A) 


Fig. 1.— Left: Phased VI OGLE-IV light curves of MACHO 176.18833.411 with the Fourier fit over-plotted. Middle: Line-of-sight radial 
velocity versus p ulsational phase for our two observations, over-plotted on a fundamental mode RRL radial velocity template and scaled by 
its l/-amplitude )|Liu|199T||. Right: Wavelength-calibrated spectrum of MACHO 176.18833.411. The CaT lines are labeled for reference. 


TABLE 1 

Properties of MACHO 176.18833.411 


a (J2000) 

<5 (J2000) 

I 

b 

<V> 

<I> 

V-amp 

Period 

Time of max brightness 
[Fe/H] 

E(F-/)min 

(m - M)o 

Heliocentric Distance 
Galactocentric Distance 
V/OS 

Vgrv 

i^a*cos(5) 

h-S 

X 

Y 
Z 

u 

V 

w 


18:00:13.08 

-27:15:39.1 

3.0795 

-1.9209 

17.586 mag 

15.969 mag 

1.27 ±0.05 mag 

0.51521996 ± 0.00000005 d 

2 456 000.21386 

—1.62±0.2 (from light curve) 

1.79±0.01 

1.21±0.03 

1.49 

14.29±0.13 mag 
7300±600 pc 
850 pc 

-372 ± 8 km s"! 

—350 ± 8 km s“^ 

8.08 ± 0.20 ± 0.40 mas yr“^ 
4.41 ± 0.20 ± 0.40 mas yr“^ 
-1.04 ± 0.56 kpc 
0.39 ± 0.03 kpc 
-0.243 ± 0.02 kpc 
-377 ± 20 km s”! 

-262 ± 26 km s”! 

-147 ± 20 km s”! 


BRAVA-RR currently has surveyed only 94 RRLs, in¬ 
cluding MACHO 176.18833.411. Is such a bulge star 
therefore really anomalous? To give us some idea if this 
RRL in fact belongs to the bulge, we integrated its orbit 
through an assumed Galactic potential which is a sum 
of the potential of a logarithmic halo, Miyamoto-Naga i 
disk, and a Hernquist bulge, as in Hawkins et al. (2015). 
Uncertainties in the orbital integrations were estimated 
by a Monte Carlo approach, where the initial conditions 
were varied to w ithin their uncertaint ies over 100 orbital 
integrations (see Hawkins et al.||2015 for details). 

The orbital integration over the past 1 Gyr is shown 
in Figure as is the distribution of the maximum dis¬ 
tance from the Galactic plane, Z^ax) and the minimum 
and maximum distance from the Galactic center for 100 
orbital draws. Although MACHO 176.18833.411 is cur¬ 
rently close to or in the bulge, its orbit clearly suggests it 
is not confined to the bulge. The escape velocity at the 
radius of the bulge is ~650 km s“^; this means that a 


star with the high velocity of MACHO 176.18833.411 can 
still be inside the bulge. However, it is clear that MA¬ 
CHO 176.18833.411 spends most of its time well outside 
the radius of the formal bulge/bar structure. 

We next address how likely it is for a bulge star to 
have such a negative radial velocity. The Galactocentic 
radial velocity, Vqrv, distribution of 229 giants toward 
the bulge at ((,6)=(4,—2) and 320 RCGs toward the bulge 
at (Z,6)=(2.4,—2.2) is shown in Figure!^ The velocities of 
the giants are taken from APOGEE ^he Apache Point 
Observatory Galactic Evolution Experiment), which is 
part of Sloan Digital Sky Survey HI, and we analyse data 
from data release 12 (DR12, Alam et al. 2015). The ve¬ 
locities of the RCGs are taken from the GIRAFFE Inner 
Bulge Survey (GIBS, Zoccali et al. 2014). The mean 
velocity and velocity dispersion of both the APOGEE 
giant sample and the GIBS RCGs are within 1 sigma of 
the mean velocity and velocity dispersion for the RRLs 
in the (^,6)=(3,—2.5) BRAVA-RR field (Kunder et al. 
2016, in prep). Assuming the radial velocity distribu¬ 
tion for the giants, red clump and RRLs is Gaussian, a 
star with VG!/jy=—350 km s“^ is a four-sigma outlier in 
velocity space. _ 

The ^9000 BRAVA giant sa mple dKunder et al. 12 0121), 
H7000 ARGOS RCG sample ( Free^n et al.||20'^ and 


the ^1200 Gaza-ESO S urvey bulge RCG sample (Rojas- 


Arriagada et al. 2014) probe further from the Galac¬ 
tic plane than where MACHO 176.18833.411 resides, so 
more contamination from e.g., the disk and halo would 
be expected, but these large radial velocity surveys can 
also give an indication of how unusual it is for a star to¬ 
ward the bulge to have such a negative velocity. Within 
the Gom-ESO Survey bulge sample, there is no star with 
Vgrv less than —350 km s“^ and <0.1% of stars have 
Vgrv less than —300 km s“^. For both the BRAVA and 
ARGOS samples, ~0.02% of stars have Vgrv less than 
—350 km s“^ and ^0.1% of stars have Vgrv less than 
—300 km s“^. 


4.2. A halo star 

There are a few possibilities that can explain the pres¬ 
ence of MAGHO 176.18833.411. The most likely is that 
this star is a halo interloper that happens to be at the 
same distance and location as where bulge RRLs re¬ 
side. This is supported by our orbital solution of MA¬ 
CHO 176.18833.411, which shows its Rmax is mostly 
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Fig. 2.— Left: A 1 Gyr orbital integration for MACHO 176.18833.411 (thick line). The open circle represents the Sun, the black ’x’ 
represents the Galactic center and the black asterisk represents the current position of the RRL. The thin grey lines are 100 draws of the 
orbital integration to illustrate the uncertainty of the orbit. MACHO 176.18833.411 has an orbit consistent with that of the halo. Right: 
The distribution of Rmim Rmax and Zmax of the 100 orbital draws for our RRL. 
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Fig. 3.— The Galactocentric velocity {Vgrv) distribution of 
APOGEE giants (solid), BRAVA giants (dashed) and GIBS RCGs 
(dashed-dot) situated spatially close to MACHO 176.18833.411 (ar¬ 
row). The mean GRV, dispersion of the Vqry and total number 
of stars are also given. 


larger than 8 kpc, strongly suggesting this RRL is not 
confined to the bulge. The time spent in the bulge is 
145±100 Myr, which corresponds to less than 15% of the 
integration. The Zmax of ~8 kpc suggests this star is not 
part of the disk, which is also consistent with its high ve¬ 
locity. This star is on an elliptical orbit (orbital ellipticity 
is 0.95 with a period of ^140 Myr) and resembles a star 
in the halo. 

Numerical simulations which match the kinematic ob¬ 
servations of the bulge well do not predict s uch high ve¬ 


locity bulge stars in this region of the sky (Shen et al. 


20101, although we note that it is possible tor signih- 


cant radial velocity outliers to be present in any popula¬ 


tion. Instead, high velocity RRLs have been found in the 
halo (e.g., ~1.5% of the local RRLs in the Layden 1994 
and Kollmeier et al. 2013 samples have Vgrv less than 
—300 km s“^, which is more than a factor of 10 larger 
than high velocity giants and RCGs found in the bulge 
surveys discussed previously). Similarly within the halo 
globular cluster (GC) sample, 3% have Vqrv less than 
—350 km s“^. 

If this RRL is a halo interloper, we presume additional 
halo contamination exists in the bulge RRL sample, with 
velocities indistinguishable from bulge stars. We would 
then expect at least a few percent of the RRLs located 
in the direction toward the bulge to be halo constituents. 
To obtain a rough approximation of the presence of the 
halo, the relation of the RRL density pr ofile in the halo 


is ex trapolated to the Galactic center (Watkins et al. 
20091 as shown in Figure To find the total number 


of halo RRLs expected in our observations, we use this 
relation along with the following approximations: (1) the 
bulge is a cylinder with a radius of 2 kpc and a minor-to- 
major axial ratio of 0.55, and (2) the bulge extends from 
I = —10 to I =4-10 and b = —10 to b =4-10 degrees, and 
so the AAOmega 3 square degree field covers 3/314 of the 
Galactic bulge. It follows, then, that ^85±70 halo RRLs 
are expected with a Galactocentric distance of 0.9 kpc in 
a 3 square degree field. The uncertainty in the number 
of halo RRLs is a function of the exact scale height and 
length of the bulge, and is illustrated in Figure The 
total number of OGLE RRLs in our 3 square degree field 
is ^1400, and the Galactocentric distance of these stars 
peaks at 0.8 kpc. Therefore, assuming the OGLE sample 
of RRLs is relatively complete, ~6% of the OGLE RRLs 
stars toward {l,b)={3,—2.5) are expected to be from the 
halo, and it is not unexpected for a few percent of halo 
stars to exist so close, both in spatial and in distance 
distribution, to the Galactic center. Proper motions of 
a larger sample of RRLs would be desirable to identify 
halo RRLs from bulge RRLs. 

It is worth noting that the orbit of MA- 
GHO 176.18833.411 has its largest excursions per¬ 
pendicular to the plane - this extreme high velocity star 
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Fig. 4. — Top: The RRL number density as a function of galac- 
tocentric distance, RqCi the 0.75 - 1.5 kpc rang e using a power 
law derived from SDSS Stripe 82 RR Lyrae stars ( [Watkins et al.| 
|2009||, where (no,Ro, ct) = (0.26 kpc“^, 23 kpc, 2.4). Middle: The 
total number of halo RRLs, NhaioRRL> within the approximate 
bulge volume as a function Rqc- The thick grey lines indicates 
the Nhaio RRL assuming a bulge radius of 2 kpc and a minor to 
major axis ratio of 0.55; the thin black lines indicate N^aio RRL 
assuming a bulge radius of 1.5 and 2.5 kpc. Bottom: Same as the 
middle panel, for for a 3 square degree bulge field. 

is confined (statistically) vertically. If the most extreme 
halo stars have an almost spheroidal zone of excursion, 
the (halo) population from where this star originated 
from is not the most extreme. It may well be that this 
RRL is a relic from an earlier era, but likely not the 
earliest. 


4.3. A bulge star 

A second possibility is that MACHO 176.18833.411 
is on the extreme edge of the bulge RRL radial veloc¬ 
ity distribution. This could be in li ne with the no¬ 
tion discussed by Dekany et al. (2013), that the bulge 
RRLs follow a spheroidal, centrally concentrated dis¬ 
tribution, as then eccentric orbits and therefore stars 
with large radial velocities, would be expected. If the 
bulge RRLs follow a Gaussian radial velocity distribu¬ 
tion with a tT=200 km s“^, it would be expected to find a 
bGi?v=—350 km s“^ radial velocity star. Such a large ve¬ 
locity distribution is not observed for the 94 BRAVA-RR 
stars in this area of the sky, but although a Gaussian dis¬ 
tribution is convenient, it is often wrong out on the tails, 
in a sense that the tails are too small. It has also been 
shown that multiple RRLs populations exist in the bulge 
( Pietrukowicz et al.||2014 |, and so there may be a small 
population ot RRLs more kinematically hot than the ma¬ 
jority of bulge RRLs, and that MAGHO 176.18833.411 
is part of this more kinematically hot component. 

Another motivation for MAGHO 176.18833.411 be¬ 
ing a member of the bulge RRL population is that a 


similar proportion of high negative velocity stars are 
found in bulge stars exhibi ting maser emissions (in both 


OH/IR and SiO masers, van Langevelde et al. 1992 
Fuji! et al.| 2006). However, it is unclear why these 


mass-losing infrared objects, typically Asymptotic Gi¬ 
ant Branch (AGB) st ars, which are thought to b e much 
younger than RRLs, ( Mouhcine fc Langon |200'2 ^1 to a 
few Gyr, ), would be more dynamically similar to bulge 
RRLs than red giants or RCGs. 


4.4. A runaway star 

A third scenario to explain the velocity of MA¬ 
CHO 176.18833.411 is that it was ejected from the Galac¬ 
tic center due to an interaction - e.g., a star-binary or 
a star - black hole collision, or ejected out of a globu¬ 
lar cluster. It is unlikely that this star is a hypervelocity 
star (HVS), as its velocity is still wel l within the realm of 
being bound to the Milky Way (e.g., Kenyon et al.|2008 l 
MACHO 176.18833.411 may be a runaway star, th ough. 


Evolved sta rs have been shown to be runaways ( Kilic 


et al. 


et al. 


2013), as have horizontal branch stars (Pereira 


2013), although by far the largest population of 


runaway stars currently known are younger, more mas¬ 
sive stars (e.g. Bromley et al.|[2009 ). Runaway stars are 
expected to dominate at low Galactic latitudes and are 
located preferentially in the di rection of the Galacti c cen¬ 
ter, between Z=325 - 35 (e.g., Bromley et al.|[2009 ). The 
expected metallicity of runaways is ^[I'e/HjL— 1.5 dex, 
slightly more metal-rich than the halo and comparable 
to MACHO 176.18833.411. 


5. CONCLUSIONS 

In this paper, we take a detailed look at a high 
negative-velocity RRL observed toward the Galactic 
bulge found in the BRAVA-RR survey. Stars of such 
high velocity are rare among bulge giants, and the more 
precise distance of the RRL makes it possible to explore 
its origin in greater detail, by integrating its orbit. We 
argue that MACHO 176.18833.411 is most likely a halo 
interloper, suggesting that contamination from halo stars 
is relevant when attempting to trace out the metal-poor 
tail of the bulge’s metallicity distribution function (e.g., 
Howes et al.||2014 ). 
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